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The  numerical,  three-dimensional,  wind-driven  circulation  study 
of  the  effect  of  a proposed  offshore  airport  on  the  thermally  stratified 
lake  conditions  in  Lake  Erie  near  Cleveland,  Ohio,  was  sponsored  by  the 
Lake  Erie  Regional  Transportation  Authority  (LERTA) , Cleveland,  Ohio, 
as  a part  of  the  model  feasibility  investigation  being  conducted  at 
the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  . The  WES 
investigation.  Task  17  of  the  LERTA  investigation,  is  a portion  of  the 
third-phase  airport  feasibility  study  undertaken  by  LERTA  to  evaluate 
airport  sites,  one  of  which  is  in  Lake  Erie  near  Cleveland. 

This  report  and  the  numerical  analyses  were  prepared  and  conducted 
by  Dr.  Donald  L.  Durham  of  the  Wave  Dynamics  Division  (WDD),  WES,  and 
Dr.  Donald  C.  Raney,  who  was  working  with  WDD  while  on  loan  from  AME 
Department,  University  of  Alabama  through  the  Intergovernmental  Person- 
nel Exchange  Program.  This  study  was  conducted  under  the  general  super- 
vision of  Dr.  R.  W.  Whalin,  Chief  of  the  Wave  Dynamics  Division,  and 
Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory.  Assisting  in 
data  reduction  for  this  study  were  Mr.  K.  A.  Turner,  Mrs.  Judy  Jones, 
and  Mr.  R.  E.  Ankeny,  WDD.  The  authors  express  their  appreciation  to 
Dr.  John  Paul,  who  is  working  with  Large  Lakes  Research  Station,  U.  S. 
Environmental  Protection  Agency  while  employed  with  Case-Western  Re- 
search University,  for  his  assistance  in  supplying  revised  computer 
codes  and  his  suggestions  for  code  conversion. 


Directors  of  WES  during  the  conduct  of  this  investigation  and  the 
preparation  and  publication  of  this  report  was  COL  G.  H.  Hilt,  CE,  and 
COL  John  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


m 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 

By 

To  Obtain 

inches 

2.54 

centimetres 

inches  per  second 

2.54 

centimetres  per  second 

2 

inches  per  second 

6.4516 

2 

centimetres  per  second 

pounds 

4.448  x 105 

dynes 

feet 

0.3048 

metres 

feet  per  second 

0.3048 

metres  per  second 

miles 

1.6093 

kilometres  per  hour 

miles  per  hour 

1.6093 

kilometres  per  hour 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

| 

* To  obtain  Celsius  (C)  temperature  reading  from  Fahrenheit  (F)  reading 
use  the  following  formula:  C ■ (5/9) (F  - 32).  To  obtain  Kelvin  (K) 
readings,  use:  X - (5/9) (F  - 32)  ♦ 273.15. 
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LAKE  ERIE  INTERNATIONAL  JETPORT  MODEL  FEASIBILITY  INVESTIGATION 


RESULTS  OF  NUMERICAL,  THREE-DIMENSIONAL,  WIND-DRIVEN  CIRCULATION 
ANALYSIS  FOR  THERMALLY  STRATIFIED  LAKE  CONDITIONS 

PART  I : INTRODUCTION 


1.  The  Lake  Erie  Regional  Transportational  Authority  (LERTA) 
is  conducting  a feasibility  and  site  selection  study  for  a major  hub 
airport  in  the  Cleveland  Service  area.  One  of  the  possible  sites  being 
evaluated  is  an  offshore  site  in  Lake  Erie  near  Cleveland,  Ohio.  As 
a part  of  the  feasibility  analysis  of  an  offshore  site,  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  conducted  a model  feasi- 
bility investigation  and  is  performing  numerical  model  studies.  The 
results  of  these  efforts  are  being  published  in  a series  of  reports 
under  the  general  title  "Lake  Erie  International  Jetport  Model  Feasi- 
bility Investigation.”  This  miscellaneous  paper  presents  preliminary 
results  of  the  numerical  simulation  analyses  of  wind-driven  circulation 
for  thermally  stratified  lake  conditions  in  Lake  Erie.  These  data  are 
results  from  WES'  converting  and  revising  a set  of  computer  programs 
originally  developed*  at  Case-Western  Reserve  University  and  WES  running 

a 12-mph  wind  from  a west  direction.  This  steady-state  wind  field 

2 

represents  the  modal  wind  speed  during  the  summer  months  (June-August) 
and  the  direction  of  the  steady-state  wind  field  possessing  the  maximum 
pertru*ation  produced  by  the  proposed  jetport  island  on  the  steady- 
state,  wind-driven  circulation^  for  well-mixed  lake  conditions 
in  Lake  Erie. 


PART  II:  HYDRODYNAMIC  MODEL 


Mathematical  Formulation 

2.  As  part  of  WES  numerical  model  feasibility  study, * the  scheme 
which  is  used  in  this  investigation  was  selected  from  a limited  number 
of  models  which  calculate  wind-driven  circulation  for  thermally  strati- 
fied lake  conditions  in  large  lakes.  The  selected  model  was  developed 
and  initially  applied*  in  testing  the  effects  of  a proposed  jetport 
island  offshore  of  Cleveland  on  the  hydrodynamics  of  the  nearshore 
region  containing  the  island.  Details  of  this  numerical  model  can  be 
found  in  References  1,  4,  and  5.  An  abbreviated  statement  of  the 
hydrodynamic  model  and  assumptions  are  presented  here  in  summary  form. 

3.  Basic  equations  for  the  numerical  model  are  derived  from  the 
time-dependent,  three-dimensional  equations  of  motion  for  a viscous, 
heat-conducting  fluid.  Figure  1 is  a schematic  of  the  model  geometry. 

In  deriving  the  model,  the  following  assumptions  are  made: 

a.  Pressure  is  assumed  to  vary  hydrostatically;  therefore, 

3P 

SI"  p*- 

b.  The  rigid-lid  approximation  is  made,  i.e.,  w(z=0)=0. 

£.  The  Boussinesq  approximation  which  assumes  that  density 
variations  are  small  and  can  be  neglected  except  in  the 
gravity  term  is  made. 

d.  Heat  sources  and/or  sinks  in  the  fluid  are  neglected. 

£.  Eddy  coefficients  are  used  to  account  for  the  turbulent  and 
molecular  diffusion  effects  in  both  the  momentum  and  energy 
equations.  The  horizontal  coefficient  is  assumed  to  be 
constant  but  the  vertical  coefficient  is  assumed  to  be 
dependent  on  the  local  vertical  temperature  gradient. 

Jf.  Variations  in  bottom  topography  are  assumed  to  be  gradual. 

4.  This  numerical  model  allows  for  variations  in  the  depth  of  the 
lake  basin.  It  uses  a nonconformal  mapping  procedure  to  stretch  the 
vertical  coordinate  with  respect  to  the  local  depth  h(x,y).  The  basic 
hydrodynamic  equations  are  transformed  according  to 
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x +-*■  X, 

y y, 

a -*-*  z/h(x,y). 

The  equations  to  be  solved  are  more  complicated  looking  because  of  the 
appearance  of  the  depth  in  the  equation,  but  they  are  solved  for  a basin 
of  constant  depth  in  the  transformed  system  which  greatly  reduces  the 
programming  complexities  of  the  model  and  makes  the  inclusion  of  depth 
variations  simpler.  A reduced  form  of  the  transformed  diffusion  terms 
are  used  by  assuming  the  terms  containing  derivatives  of  the  depth  are 
neglected  with  respect  to  those  terms  containing  only  the  depth.  This 
approximation  is  used  in  meteorological  problems  when  topographic  vari- 
ations are  included  (Refs.  6 and  7).  The  bottom  topography  h(x,y) 
in  the  nearshore  region  of  Lake  Erie  around  Cleveland,  Ohio,  were 
obtained  from  Lake  Survey  Center  charts  of  Lake  Erie.  An  outline  of  the 
procedure  with  which  these  values  were  read  from  charts,  interpolated 
and  smoothed  is  given  in  Reference  8. 

5.  The  resulting  system  of  transformed  equations,  in  non-dimensional 
form,  are  the  following: 


i 3(hu)  t i.  a(hv)  t an  _ 0 

h »x  h 3y  do 


S * -I  S ^ ♦ i - & * i «•&  ♦ 

- “»>'  * £>’  k <*g>. 

- ♦ Re | » iCSHT)  ♦ i »Si)  _ Rou  - - f?  ♦ 1 ♦ A-  <„»_*), 

*J  h h h iy  3o  1 »y  h *3x  »x'  ay  'nay'* 

- ^ ih/ifdo + U (.fa^do  - aap)i ♦ <&>’  i ^ 
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\ 
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M*AT+  R.l  a(huAT)  1 a(hvaT)  + anaj,  M i ,2.  (haAT.  * _a_  fhaAl>t 

' >L  1.  k lu  1_  1 k 1 i..  ' i..  ' 'I 


ax 


where: 


w = 


u 

V 

Wbn 

u0h0 


b0* 


r - PRe 

p^ghoFr7’ 

at=— Te 
te  • 


7 = 


= *¥ 
ah’ 

; = *b0j 


Ro  = 

ah  ’ 

Pr  = ^H 

BH’ 

Po  = P(TE)  and 


h ay 

(r*) 


do 

-f  f^*)1  _L  2 /»®at. 


h 8x  8y  ay  dy 


h. 


n 


h*  a^’ 


1 3P  _ Re  ..  . 

rr  77  '*  + *<»)• 

h 3o  Fr 


Ap  = f(AT). 
o - z/h(x,y). 


I |w  - o(u— - + v^)l  = 42, 

h ax  ay  1 dt 


U0 

- JL 

b0’ 


ho’ 

Ap  = P ~ Pa. 
_ 

. = *ah 
1?’ 


0 = 


By 

BH’ 


Fr  = — VP- 

s/fro  ’ 

Re  = 

AH 


uq  = reference  velocity, 

b = horizontal  reference  length,  . 
o 

hQ  = vertical  reference  length, 

Ajj  = horizontal  eddy  viscosity. 

Ay  = vertical  eddy  viscosity, 

Bu  = horizontal  eddy  diffusivity. 

By  = vertical  eddy  diffusivity, 

T_  « equilibrium  temperature  which  is  temperature  at  the  surface 
for  which  there  is  no  heat  transfer. 

k = Coriolis  parameter, 

f (AT)  = equation  of  state  and 

(”)  ■ refers  to  dimensional  quantity. 
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6.  The  Poisson  equation  for  pressure,  which  contains  the  rigid- 
lid  condition,  is  derived  by  taking  the  divergence  of  the  vertically 
integrated  horizontal  momentum  equations  and  using  the  vertically 
integrated  continuity  and  hydrostatic  pressure  equations.  The  Poisson 
equation  is 


. 3P  a 3P 

h h 33T + h h = -h  It 

l)  2 

, o.  3 ,1  3u.  1 3u. 

* (h->  ta  (h  Y - h * > 

o 0=1  a=0 


,bo,2  3 ,1  3v,  1 3v.  . 

(h  } 3y  (h  Y 3o5  " h Y 3o5  5 


o=l 


o=0 


3 L A m A 3hu  . 1 3huv  3Jlu^  „ 1 3 t3u 

- * h ' [Re(h  IT  * iT  aj—  * iT>  * Rov  - h 57 


1 3 v,  . Re  ( 3 i r A , 3h  A v , 

h 3?  h 37  72  (37  h 1 Apd°  - ° 37  Ap)]d° 

Fr  o 


3 t.A  rr,.A  3huv  . 1 3hv  . 3ftv.  ' 13  ^ 3v 

37  H [Re(h  -IT  + h IT"  + 37"J  ' R0U  ' h 37  h 37 


JL  — v,  — 
h 3y  3y 


Fr 


r-  r-  h rrr  * (§7  h£  Apdo  " 0 f~  Ap)]do 


where  Pg  is  the  integration  constant  resulting  from  the  vertical  inte- 
gration of  the  hydrostatic  pressure  equation  and  is  the  surface  pressure 
i.e.,  the  pressure  at  the  surface  z = 0. 

7.  The  following  boundary  conditions  are  used  with  the  above  system 
of  equations. 

u = gjfy.z) 

River  outflow  v = gjCy.z) 

AT  = g3(y,z) 
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Shore 


u = 0 
v = 0 


Bottom 


u = 0 
v = 0 
w = 0 


Surface 


= twx 


3v 

Km  twy 


= KAT 


Other  Boundary 
either 


w = 0 


u = f. 


0 

3y 


v = f„ 


Pressure  Conditions 


AT  = f. 


— = integrated  x or  y momentum  equation, 
specify  pressure  level  at  one  point. 


The  functional  forms  g^,  g£,  and  g^  are  the  specified  velocity  and 
temperature  profiles  across  the  river  outfall.  Boundary  conditions 
at  the  outer  x and  y boundaries  are  either  that  the  normal  derivatives 
of  the  velocity  and  temperature  are  zero,  or  that  the  velocity  and 
temperature  are  specified  (fj,  fj,  fj) . 


10 


Numerical  Procedure 


8.  The  general  arrangement  of  variables  in  the  numerical  grid 
system  is  identical  to  that  used  previously  by  Paul  and  Lick1.  Horizon- 
tal velocities  are  defined  at  integral  nodal  points,  temperature  is  de- 
fined at  half-integral  nodal  points  in  the  horizontal  and  integral  nodal 
points  in  the  vertical,  and  the  surface  pressure  is  defined  at  half- 
integral nodal  points  in  the  horizontal.  Figures  2 and  3 are  the 
horizontal  grid  and  typical  vertical  grid  sections  for  the  nearshore 
model.  The  relative  positions  of  the  various  variables  within  the 
numerical  grid  are  depicted  in  these  figures.  The  location  of  the 
jetport  is  indicated  by  shaded  cells  in  Figure  2. 

9.  The  finite  difference  approximations  to  the  equations  are 
derived  by  integrating  the  equations  over  nodal  cells  (Figure  4)  using 
either  the  mid-point  or  trapezoidal  integration  rule  to  evaluate  these 
integrals.  In  the  derivation  of  the  finite  difference  equations,  a 
simple  average  of  neighboring  values  is  used  for  variables  which  are 
not  defined  at  required  points.  The  Euler  explicit  time  scheme  is  used 
exclusively  in  the  present  model.  Details  of  the  finite  difference 
approximation  to  the  hydrodynamic  equations  are  given  in  References  1 
and  9. 

10.  Solution  of  the  difference  equations  is  obtained  by  the 
following  scheme: 

a.  Values  from  the  previous  time  step  are  assumed  to  be 
available. 

Ik  Temperature  is  calculated  by  an  explicit  time  scheme. 

£.  Density  is  calculated  from  the  equation  of  state. 

d.  Surface  pressure  is  calculated  with  the  right -hand -side 
of  the  equation  evaluated  from  the  new  temperatures  and 
density  values  and  other  previous  time  step  values. 

e_.  Horizontal  velocities  are  calculated  by  an  explicit  time 
scheme. 

f.  Vertical  velocity  is  calculated  by  vertically  integrating 
the  continuity  equation  from  the  bottom. 

£.  The  present  time  step  is  now  complete. 
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Figure  2.  Location  of  variable*  within  the  horizontal  grid  for  the 
nearshore  model.  Jet port  island  1*  indicated  by  shaded  cells 
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Figure  3<  Location  of  variables  in  vertical  sections  for 
the  nearshore  model 
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11.  The  temperatures  are  calculated  by  the  explicit  time  scheme 

and  are  checked  for  static  stability,  i.e.,  if  temperatures  decrease 

monotonically  downward  (assuming  that  density  increases  with  decreasing 

temperature).  When  static  instabilities  are  encountered,  an  infinite 
. . 9 

mixing  procedure  is  used  with  temperatures  over  any  unstable  region 
being  averaged.  Using  temperature  values  at  new  time  steps,  density  is 
calculated  from  the  equation  of  state  of  fresh  water  which  is  assumed 
to  vary  linearly  with  temperature. 

12.  At  each  time  step,  the  Poisson  equation  for  the  surface  pres- 
sure is  solved  by  the  point  successive-over-relaxation  (SOR)  method.  In 
the  Poisson  equation,  the  forcing  term  involves  a time  derivative  of  the 
vertical  velocity  at  the  surface.  The  vertical  velocity  at  the  surface 
is  zero  by  the  rigid-lid  condition.  However,  non-zero  values  for  the 
vertical  velocity  are  obtained  numerically  by  vertical  integration  of 
the  continuity  equation  and  indicate  that  the  continuity  equation  can 
not  be  satisfied  exactly  by  the  finite  difference  solution.  This  error 
can  grow  in  time;  thus,  the  Hirt-Harlow  corrective  procedure1  is  used 
to  correct  for  this  error. 

Nearshore  Application 

13.  A particular  finite  difference  grid  depends  on  the  actual  ge- 
ometry to  be  described.  The  shoreline  used  in  this  model  is  determined 
by  the  20-ft  contour.  Zero  depth  is  not  chosen  as  the  shoreline  because 
the  vertical  coordinate  transformation  is  singular  for  zero  depths. 

This  use  of  non-zero  depth  shorelines  does  not  appreciably  affect  the 
boundary  and  is  not  a restriction  on  the  model.  The  nearshore  region, 
modeled  in  this  study,  is  a 16-mile  by  16-mile  area  in  Lake  Erie  near 

Cleveland,  Ohio,  and  is  similar  to  the  area  used  in  previous  WES 
3 8 

studies.  * A constant  horizontal  spacing  of  one  mile  is  used  with  the 
proposed  jetport  represented  by  a two-mile  by  three-mile  island  five 
miles  off  Cleveland  in  approximately  50  ft  of  water  (Figure  2). 
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14.  The  vertical  eddy  coefficient  is  taken  as  dependent  on  the 

local  vertical  temperature  gradient.  This  form  is  similar  to  that 

11  12 

suggested  by  Sundaram,  et  al  (1969,  1970)  ’ and  is  identical  to  that 

used  in  a previous  application  of  this  model. ^ The  expression  for  the 
vertical  eddy  coefficient  Ay  is: 

Ay  = a ♦ S £ 

where  a and  $ are  constants  dependent  on  the  local  conditions  of  the 
physical  system  modeled.  The  constant  a is  chosen  so  that  in  the  ab- 
sence of  vertical  temperature  gradients,  the  eddy  coefficient  is  equal 
to  that  which  would  be  used  for  a constant  eddy  coefficient. 

15.  The  wind  stress  imposed  on  the  water  surface  due  to  the  wind 
action  is  calculated  from  the  formulae  developed  by  Wilson  (1960). 

3 4 8 

These  formulae  have  been  successfully  used  in  numerical  calculations  ’ ’ 
of  wind-driven  circulations  in  lakes  and  in  a previous  application1"*  of 
this  present  model  to  a power  plant  outfall. 

16.  Boundary  conditions  for  the  open  water  boundaries  of  the 
model  are  as  follows: 

a.  Along  the  outer  x boundary  in  the  lake  (boundary  3 of  Figure  5), 
velocities  and  temperature  are  specified.  Pressure  is  obtained  from 

the  vertically  integrated  x momentum  equation.  These  boundary  values 
of  velocity  and  temperature  were  inferred  using  good  engineering  judg- 
ment from  an  application1  of  the  numerical  model  to  either  the  Central 
Basin  of  Lake  Erie  or  the  entire  lake  and  WES  results^  of  the  steady- 
state,  wind-driVen  circulation  study. 

b.  Along  the  two  y boundaries  in  the  lake  (boundaries  2 and  4 of 
Figure  5),  the  variables  are  assumed  to  be  smoothly  varying,  i.e.,  the 
first  normal  derivatives  of  the  velocities  and  temperature  are  zero 
and  the  second  normal  derivative  of  the  pressure  is  zero. 

17.  In  the  non-dimensional  formulation  of  the  governing  equa- 
14  5 

tions,  ’ * various  characteristic  parameters  are  defined  and  are  used 
in  scaling  the  non-dimensional  results  of  the  numerical  study.  A list 
of  the  pertinent  parameters  and  their  values,  used  in  this  study,  is 
presented  in  Table  1. 


J 
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BV  = \ 

BH  = *H 

Uo 

k 

epilimnion  temperature 
hypolimnion  temperature 
thermocline  depth 
wind 


Table  1 

in  Application  of  Model 


16  mi  (2.58  x 106  cm) 

80  ft  (2440  cm) 

aT  ? 

[16.8  + (5.0)  cm  /sec 
5 x 103  cm^/sec 

15  cm/sec 

in-4  -1 
10  sec 

75°f  (24°C) 

55°f  (13°C) 

30  ft  (915  cm) 

12  mph  (536  cm/sec)  from  west 
77.4 

1.33  x 103 
9.6  x 10"3 
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1 4 

18.  The  numerical  model,  ’ which  was  developed  at  Case  Western- 
Reserve  University  to  compute  three-dimensional,  time-dependent,  wind- 
driven  currents  for  a thermally  convective  fluid  in  a nearshore  region 
of  Lake  Erie,  consisted  of  a Fortran  code  to  compute  the  velocity 
components  (u,  v,  and  w) , temperature,  and  density  at  any  selected 
depth  in  the  water  column  and  pressure  at  the  surface.  The  Waterways 
Experiment  Station  converted  this  code  to  operate  on  a CDC  7600  Computer 
and  enhanced  the  output  procedures  to  include  various  graphics  routines. 


PART  III:  RESULTS  AND  CONCLUSIONS 


Results  of  Model  Application 

19.  The  numerical  hydrodynamic  model,  described  in  Part  II,  has 
been  used  to  investigate  the  effect  of  a proposed  jetport  island  on  the 
summer  stratification  pattern  in  the  nearshore  area  offshore  of  Cleve- 
land, Ohio.  Due  to  monetary  and  time  constraints  in  the  feasibility 

study,  only  two  steady-state  wind  fields  were  considered  in  studying  the 

2 

summer  (stratified)  lake  conditions.  The  modal  wind  speed  (12  mph)  and 

2 

most  frequently  occuring  direction  (SSE-SW)  for  the  months  of  June, 
July,  and  August  were  considered  as  the  typical  steady-state  summer  wind 
field.  This  wind  direction  has  an  annual  frequency  of  occurrence  of 

42  percent  and  a trimonth  (June-August)  frequency  of  occurence  of 

43  percent.  A previous  WES  study1  contracted  to  Case-Western  Reserve 
University  used  the  above  wind  field  in  studying  the  effects  of  a pro- 
posed jetport  island  on  the  summer  circulation  near  Cleveland,  Ohio. 
Results  of  this  study  are  presented  in  detail  in  Reference  1.  For  this 
south  wind  direction,  which  produced  minimum  effects  in  the  steady- 
state,  wind-driven  circulation  for  well-mixed  (fall-winter)  lake  condi- 

3 

tions  , the  effects  of  the  proposed  jetport  island  on  the  velocity  re- 
gime were  localized  within  2 to  3 miles  of  the  island.  For  the  temper- 
ature regime,  the  effects  of  the  proposed  jetport  island  extend  several 
miles  from  the  island;  however,  these  effects  did  not  extend  to  the 
shoreline  near  Cleveland  due  to  the  particular  direction  (south)  of  the 
wind  field.  In  addition  to  the  above  application  (12  mph  south  wind), 
WES  chose  another  steady-state  wind  field  in  studying  the  summer  lake 
conditions.  A 12  mph  west  wind  was  chosen  for  the  second  wind  field 

3 

since  this  direction  produced  maximum  effects  in  the  WES  study  of 
steady-state,  wind-driven  circulation  for  well-mixed  (fall-winter)  lake 
conditions.  Choosing  a 12-mph  modal  wind  speed  and  two  wind  directions 
(S5W)  which  produced  the  minimum  and  maximum  jetport  island  effects  in 
a previous  WES  study^  of  well-mixed  lake  conditions,  the  extreme  effects 
of  a jetport  island  on  wind-driven  circulation  for  stratified  lake 
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conditions  should  be  estimated  sufficiently  for  this  engineering  feasi- 
bility study. 

20.  Results  for  the  nearshore  application  of  the  hydrodynamic 
model  without  the  jetport  island  for  a 12  mph  west  wind  are  shown  in 
Plates  13  to  24,  and  results  with  the  jetport  island  are  shown  in 
Plates  25  to  36.  In  addition,  vector  plots  of  differences  in  horizontal 
velocity  regimes  with  and  without  a jetport,  contours  of  differences  in 
velocity  magnitude  with  and  without  a jetport,  and  contours  of  differ- 
ences in  temperature  regimes  with  and  without  a jetport  are  presented  in 
Plates  37  to  42,  Plates  43  to  48,  and  Plates  49  to  54,  respectively. 

The  difference  in  horizontal  velocity  and  temperature  are  computed  by 
subtracting  at  each  grid  point  in  the  nearshore  region  the  velocity  or 
temperature  with  a jetport  from  the  velocity  or  temperature  without  a 
jetport.  For  contours  of  differences  in  velocity  magnitude,  the  abso- 
lute values  of  the  magnitude  differences  are  contoured.  For  contours  of 
temperature  differences,  a positive  value  indicates  a higher  temperature 
without  the  jetport  and  a negative  value  indicates  a lower  temperature 
without  the  presence  of  a jetport.  Results  are  presented  after  7.4  hrs 
of  real  time  simulation.  The  initial  conditions  of  horizontal  velocity 
and  temperature  fields  for  the  12  mph  west  wind  are  indicated  in 
Plates  1 to  6 and  7 to  12,  respectively. 

21.  A comparison  of  horizontal  velocity  plots  (Plates  13-20,  25- 
30,  37-42)  indicate  the  proposed  jetport's  major  effects  on  the  horizon- 
tal velocity  regime  for  a steady-state  12  mph  west  wind  are  localized 
within  an  area  3 to  4 miles  from  the  jetport  island.  The  areal  extent 
and  magnitude  of  change  in  the  horizontal  velocity  regime  that  are  as- 
sociated with  the  proposed  jetport  are  vividly  depicted  in  Plates  37  to 
48.  In  Plates  43  to  48,  jetport  perturbation  of  0.25  ft/sec  in  the 
horizontal  velocity  regime  extend  to  shore  in  the  immediate  vicinity  of 
Cleveland,  Ohio.  These  0.25  ft/ sec  perturbations,  which  extend  to  shore, 
occur  throughout  the  water  column.  The  0.25  ft/sec  perturbations  in 
horizontal  velocity  extend  as  much  as  8 miles  lakeward  of  the  jetport. 
Velocity  changes  of  0.5  ft/sec  extend  a maximum  of  4 miles  from  the  jet- 
port island  and  come  within  one  mile  of  the  shoreline  near  Cleveland, 
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Ohio.  Velocity  changes  of  1 ft/sec  extend  a maximum  of  two  mile*  from 
the  jetport  island  with  the  maximum  change  of  1.5  ft/sec  occurring  at 
20  ft  depth  in  water  column.  The  major  influence  of  the  proposed  jet- 
port  island  on  the  temperature  structure  in  the  lake  extends  5 to  6 
miles  from  the  jetport  island.  The  magnitude  and  areal  extent  of  the 
jetport's  effect  on  the  horizontal  isotherms  are  depicted  in  Plates  49 
to  54.  In  addition,  these  contours  of  temperature  differences  vividly 
depict  areas  of  upwelling  of  cool  water  and  downwelling  of  warm  water 
around  the  jetport.  Temperature  changes  of  0.5°  (F)  occur  along  the 
shoreline  in  the  immediate  vicinity  of  Cleveland,  Ohio.  Temperature 
changes  ot  1.0°  (F)  extend  within  a half  mile  of  the  shoreline  and 
5-1/2  miles  lakeward  of  the  jetport  island.  Temperature  changes  of 
5.0°  (F)  or  greater  are  confined  within  2-1/2  miles  from  the  island. 

The  largest  decrease  (14°F)  in  temperature  occurs  at  the  water  surface 
in  the  strong  upwelling  region  along  the  southern  side  of  the  jetport 
island.  Large  increases  (16°F)  in  temperature  occur  at  40-ft  depth  in 
the  water  column  in  the  strong  downwelling  region  along  the  northern 
side  of  the  jetport  island. 

22.  The  large  area  of  influence  by  the  proposed  jetport  island  in 
temperature  regime  and  horizontal  velocity  can  be  associated  with  the 
upwelling  and  downwelling  regions  which  change  the  stratification  struc- 
ture in  the  lake  around  the  jetport.  Being  a variable-density  model, 
changes  in  the  temperature  structure  cause  changes  in  the  velocity 
patterns.  Therefore,  these  strong  upwelling  and  downwelling  regions 
produce  mixing  between  the  hypolimnion  and  the  epilimnion  in  this  region 
of  the  lake  and  can  result  in  the  erosion  of  the  thermocline  in  this 
area.  The  results  of  this  model  are  dependent  of  the  wind  direction 
with  shape,  size,  and  location  of  upwelling  and  downwelling  depending  on 
wind  speed  and  direction.  In  addition,  the  area  of  influence  by  the  jet- 
port would  be  expected  for  a constant  wind  field  to  change  with  time  (in- 
creasing in  areal  extent  into  lake  with  time).  Values  (Table  1)  for 

various  parameters,  in  particular  the  vertical  eddy  diffusivity  coeffi- 

14 

cient  , have  a large  effect  on  the  vertical  temperature  structure  which, 
in  turn,  influences  the  velocity  regime.  Further  experimentation  with 


22 


the  numerical  model  would  be  required  to  determine  the  influence  of 
various  parameters.  Additional  confidence  in  the  model  results  could 
be  obtained  by  comparing  results  of  model  without  a jetport  with  pro- 
totype data.  However,  no  prototype  data  are  available  at  this  time 
for  such  model  verification.  Thus,  the  results  presented  in  this  report 
should  be  considered  preliminary  and  capable  of  defining  for  engineering 
feasibility  purposes  the  qualitative  effects  of  a proposed  jetport 
island  on  the  thermal  structure  and  horizontal  velocity  in  the  lake  for 
a steady-state  12  mph  west  wind. 

Conclusions 


23.  General  conclusions  of  WES  numerical  model  study  of  the 
effects  of  a proposed  jetport  island  on  the  wind-driven  circulation  in 
Lake  Erie  for  thermally  stratified  lake  conditions  with  a 12  mph  west 
wind  are  presented  below: 

a.  The  jetport  island's  effects  of  an  engineering  and  prac- 
tical interest  are  confined  within  a 16  mile  by  16  mile 
nearshore  region  with  major  effects  in  velocity  and 
thermal  regimes  being  localized  near  the  jetport  island. 

b.  Due  to  the  presence  of  a jetport  island,  horizontal  vel- 
ocity changes  of  0.25  ft/sec  extend  as  much  as  8 miles 
lakeward  of  the  jetport  and  reach  the  shoreline  in  the 
immediate  vicinity  of  Cleveland,  Ohio.  Velocity  changes 
of  0.5  ft/ sec  extend  lakeward  a maximum  of  4 miles  from 
the  jetport  and  reach  within  1 mile  of  the  shore  near 
Cleveland.  Velocity  changes  of  1 ft/sec  or  greater  are 
confined  within  2 miles  of  the  jetport.  The  maximum 
change  in  horizontal  velocity  is  estimated  to  be  1.5  ft/ 
sec. 

£.  The  jetport' s effects  on  thermal  structure  extends  lake- 
ward  5 to  10  miles  from  the  jetport  with  0.5°  (F)  changes 
in  temperature  reaching  the  shoreline  near  Cleveland, 

Ohio.  Temperature  changes  of  1°  (F)  extend  within  a half 
mile  of  the  Cleveland  shoreline  and  5 to  6 miles  lakeward 
of  the  jetport.  Changes  of  5°  (F)  or  greater  are  confined 
within  2-1/2  miles  of  the  jetport.  Maximum  temperature 
changes  of  approximately  15°  (F)  are  estimated  in  the 
areas  of  strong  upwelling  and  downwelling  along  the  jet- 
port island's  perimeter. 


d.  The  location,  size,  and  shape  of  areas  of  upwelling  and 
downwelling  around  the  jetport  island  are  dependent  on 
wind  speed  and  direction.  These  areas  induce  mixing  be- 
tween the  epilimnion  and  hypolimnion  layers,  and  can  cause 
erosion  of  the  thermocline  in  the  immediate  vicinity  of 
the  jetport  island. 

e_.  For  the  12  mph  west  wind,  the  downwelling  along  the  shore 
northeast  of  Cleveland  is  basically  not  affected. 

24.  Based  on  above  results  of  this  study,  a 2 mile  by  3 mile  jet- 
port island  located  in  Lake  Erie  at  least  4 miles  offshore  of  Cleveland, 
Ohio  will  change  the  horizontal  velocity  regime  and  thermal  structure 
in  the  immediate  vicinity  of  Cleveland,  Ohio.  Although  these  effects 
appear  to  be  small  (0.25  ft/sec  and  0.5°F),  their  impact  on  the  circu- 
lation within  the  Cleveland  Harbor  need  to  be  studied.  Effects  of  an 
engineering  interest  occur  within  3 to  5 miles  of  the  proposed  jetport 
island.  These  results  are  preliminary,  have  not  been  verified  using 
prototype  data,  and  should  be  used  in  the  feasibility  study  to  indicate 
the  qualitative  effects  of  the  proposed  jetport  island  on  the  wind- 
driven  circulation  in  Lake  Erie  for  thermally  stratified  lake  conditions 
with  a 12  mph  west  wind. 
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INITIQL  CONDITIONS  WITHOUT  JETPORT 


HORIZONTAL  VELOCITY  AT  15  FT  DEPTH 
12  MPH  WEST  WIND 

INITIAL  CONDITIONS  WITHOUT  JETPORT 


PLATE  3 


HORIZONTAL  VELOCITY  AT  HO  FT  DEPTH 
12  MPH  WEST  WIND 
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PLATE  13 
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APPENDIX  A:  NOTATION 


A^  Horizontal  eddy  viscosity 

Av  Vertical  eddy  viscosity 

bQ  Horizontal  reference  length 

B^  Horizontal  eddy  diffusivity 

By  Vertical  eddy  diffusivity 

fj  Outer  boundary  condition  for  u 

f.  Outer  boundary  condition  for  v 

fj  Outer  boundary  condition  for  AT 

f(AT)  Equation  of  state 

Fr  Froude  number 

g Gravitational  acceleration 

gj  River  outflow  boundary  condition  for  u 

g^  River  outflow  boundary  condition  for  v 

g^  River  outflow  boundary  condition  for  AT 

h Bottom  depth 

hQ  Reference  depth 

k Dimensional  Coriolis  parameter 

K Surface  heat  transfer  coefficient 

P Pressure 

Pg  Surface  pressure 

Pr  Turbulent  Prandtl  number 

Re  Reynolds  number 

Ro  Nondimensional  Coriolis  parameter 
t Time 


A1 


twx 


Surface  wind  stress  in  x direction 


twy  Surface  wind  stress  in  y direction 

T Temperature 

Tp  Equilibrium  temperature 

u Velocity  in  x direction 

uq  Reference  velocity 

v Velocity  in  y direction 

w Velocity  in  z direction 

a Constant  in  variable  Ay  term 

8 Ratio  of  vertical  to  horizontal  eddy  diffusivities  or  constant 
in  variable  Ay  term 

y Ratio  of  vertical  to  horizontal  eddy  viscosities 
AT  Temperature  difference 

Ap  Density  difference 

p Density 

pQ  Reference  density 

o Transformed  vertical  coordinate 

fi  Velocity  in  o direction 

( ) Refers  to  dimensional  quantity 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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driven  circulation  analysis  for  thermally  stratified  lake  condi- 
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from  National  Technical  Information  Service,  1978. 
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